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The scaling of the yields of heavy projectile residues from the reactions of 25 MeV/nucleon 86 Kr 
projectiles with 124 Sn, 112 Sn and 64 Ni, 58 Ni targets is studied. Isotopically resolved yield distribu- 
tions of projectile fragments in the range Z=10-36 from these reaction pairs were measured with 
the MARS recoil separator in the angular range 2.7°-5.4°. For these deep inelastic collisions, the 
velocities of the residues, monotonically decreasing with Z down to Z~26-28, are employed to char- 
acterize the excitation energy. The ratios R,2i(N,Z) of the yields of a given fragment (N,Z) from each 
pair of systems are found to exhibit isotopic scaling (isoscaling), namely, an exponential dependence 
on the fragment atomic number Z and neutron number N. The isoscaling is found to occur in the 
residue Z range corresponding to the maximum observed excitation energies. The corresponding 
isoscaling parameters are a=0.43 and /3=-0.50 for the Kr+Sn system and a=0.27 and /3=-0.34 for 
the Kr+Ni system. For the Kr+Sn system, for which the experimental angular acceptance range lies 
inside the grazing angle, isoscaling was found to occur for Z<26 and N<34. For heavier fragments 
from Kr+Sn, the parameters vary monotonically, a decreasing with Z and /3 increasing with N. 
This variation is found to be related to the evolution towards isospin equilibration and, as such, 
it can serve as a tracer of the N/Z equilibration process. The present heavy-residue data extend 
the observation of isotopic scaling from the intermediate mass fragment region to the heavy-residue 
region. Interestingly, such high-resolution mass spectrometric data can provide important infor- 
mation on the role of isospin and isospin equilibration in peripheral and mid-peripheral collisions, 
complementary to that accessible from modern large-acceptance multidetector devices. 

PACS numbers: 25.70.Mn,25.70.Lm,25.70.Pq 



I. INTRODUCTION 



The isotopic composition of nuclear reaction products 
contains important information on the role of the isospin 
on the reaction dynamics 0]. Recently, increased inter- 
est in the N/Z degree of freedom and its equilibration 
HH H HL as well as in the isospin asymmetry depen- 
dent terms of the nuclear equation of state ja, Q has 
motivated detailed measurements of the isotopic distri- 
butions of fragments with Z>2 

011 113 E2 Isotopi- 
cally resolved data in the region Z=2-8 have revealed 
systematic trends, but their N/Z properties are unavoid- 
ably affected by the decay of the excited primary frag- 
ments. It has recently been shown E2 that isospin effects 
can be investigated by comparing the yields of fragments 
from two similar reactions that differ only in the isospin 
asymmetry. In this case, the effect of sequential decay 
of primary fragments can be bypassed to a large extent. 
It has been revealed that for statistical fragment produc- 
tion mechanism(s), if two reactions occurring at the same 
temperature have different isospin asymmetry, the ratio 
R,2i(N,Z) of the yields of a given fragment (N,Z) obtained 
from the two reactions 2 and 1 exhibits an exponential 
dependence on N and Z of the form: 



*On leave of absence from the Institute of Physics of the Slovak 
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R 2 i(N, Z) = Cexp{aN + (3Z) 

where a and (5 are the scaling parameters and C 
is an overall normalization constant. This scaling be- 
havior is called isotopic scaling or, isoscaling, [lj] and 
has been observed in a variety of reactions under the 
conditions of statistical emission and equal temperature 
[3 El 111 El- (If the temperatures are different, a 
generalized isoscaling can be obtained with appropriate 
temperature corrections |17|.) The isoscaling parame- 
ters a and (3 are shown to contain information on the 
early stages of fragment formation (before sequential de- 
excitation). It remains a matter of further investigation 
to what extent they may provide information on the type 
and the details of the reaction mechanism(s). 

Up to the present, the isoscaling phenomenon has 
been systematically investigated with isotopically re- 
solved fragments not heavier than Z=8. It would be in- 
teresting to investigate the behavior of heavier fragments 
up to the region of heavy residues. Heavy residues are 
known to comprise a large fraction of the reaction cross 
section at the intermediate energy regime [l8|. However, 
efficient collection and complete characterization of the 
residues in terms of their atomic number Z, mass num- 
ber A, ionic charge state q and velocity requires the use 
of a magnetic spectrometer. At this point, it should be 
noted that the scaling behavior of heavy fragment yields 
from light particle induced reactions at relativistic ener- 
gies on separated Sn targets has been recently reported 
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|19| . In [l9j , the fragment yields were measured with ra- 
diochemical techniques. With mass spectrometric tech- 
niques, however, apart from isotopically resolved yields, 
the velocities of the fragments can be obtained with high 
resolution and can provide information on the excitation 
energy of the primary fragments. 

In the present work, isotopic scaling is investigated for 
fragments obtained in high-resolution mass spectromet- 
ric studies of projectile fragments from the reactions of 
25 MeV/nucleon 86 Kr with 124 < 112 Sn and 64 < 58 Ni targets. 
Yield ratios of fragments with Z=10-36 were obtained 
and isotopic scaling was investigated. Using the veloci- 
ties of the residues, the excitation energy of the hot pri- 
mary fragments was characterized. Finally, the values 
of the isoscaling parameters were used to provide isospin 
information on the hot emitting sources. The paper is or- 
ganized as follows. In Section II, a brief description of the 
experimental apparatus, the measurements and the data 
analysis is given. In Section III, after examiming the av- 
erage velocity and N/Z characteristics of the fragments, 
the isotopic scaling of the fragment yields is studied in 
detail. Finally, conclusions from the present study are 
summarized in Section IV. 



II. EXPERIMENTAL METHODS AND DATA 
ANALYSIS 

The present study was performed at the Cyclotron In- 
stitute of Texas A&M University. A 25 MeV/nucleon 
86 Kr 22+ beam from the K500 superconducting cyclotron, 
with a typical current of ~1 pnA, interacted with isotopi- 
cally enriched targets of 124 Sn, 112 Sn (2 mg/cm 2 ) and 
64 Ni, 58 Ni (4 mg/cm 2 ). The reaction products were an- 
alyzed with the MARS spectrometer [2(| offering an an- 
gular acceptance of 9 msr and momentum acceptance of 
4%. The primary beam struck the target at 4.0° rel- 
ative to the optical axis of the spectrometer. The di- 
rect beam was collected inside the target chamber on a 
square Faraday cup lying outside of the angular accep- 
tance of the spectrometer. Fragments were accepted in 
the polar angular range 2.7°-5.4°. This angular range 
lies inside the grazing angle of 6.5° of the Kr+Sn reac- 
tions and mostly outside the grazing angle of 3.5° of the 
Kr+Ni reactions at 25 MeV/nucleon |2l|. It should be 
noted that the spectrometer angle setting was chosen to 
optimize the production of very neutron-rich fragments 
from the Kr+Sn systems whose detailed study is reported 
elsewhere jsisj. An Al foil (1 mg/cm 2 ) was used to reset 
to equilibrium the ionic c harg e states of the projectile 
fragments. MARS optics |20j provides one intermedi- 
ate dispersive image and a final achromatic image (focal 
plane). At the focal plane, the fragments were collected 
in a 5x5 cm two-element (AE, E) Si detector telescope. 
The AE detector was a position-sensitive Si strip detec- 
tor of 110 /xm thickness, whereas the E detector was a 
single-element Si detector of 950 /xm, respectively. Time 
of flight was measured between two PPACs (parallel plate 



avalanche counters) positioned at the dispersive image 
and at the focal plane, respectively, and separated by 
a distance of 13.2 m. The PPAC at the dispersive im- 
age was also X-Y position sensitive and used to record 
the position of the fragments. The horizontal position, 
along with NMR measurements of the field of the MARS 
first dipole, was used to determine the magnetic rigidity, 
Bp, of the particles. The reaction products were charac- 
terized by an event-by-event measurement of energy-loss, 
residual energy, time of flight, and magnetic rigidity. The 
response of the spectrometer/detector system to ions of 
known atomic number Z, mass number A, ionic charge 
q and velocity was calibrated using a low intensity 86 Kr 
primary beam and other beams at 25 MeV/nucleon. 

The determination of the atomic number Z was based 
on the energy loss of the particles in the AE detector 
[23| and their velocity. The ionic charge q of the parti- 
cles entering MARS was obtained from the total energy 
E tot =AE + E, the velocity and the magnetic rigidity. 
The measurements of Z and q had resolutions of 0.5 and 
0.4 units (FWHM), respectively. Since the ionic charge 
must be an integer, we assigned integer values of q for 
each event by putting windows (Aq = 0.4) on each peak 
of the q spectrum. Using the magnetic rigidity and ve- 
locity measurement, the mass-to-charge A/q ratio of each 
ion was obtained with a resolution of 0.3%. Combining 
the q determination with the A/q measurement, the mass 
A was obtained as: A — qi nt x A/q ((lint is the integer 
ionic charge) with a resolution (FWHM) of 0.6 A units. 
Combination and normalization of the data at the var- 
ious magnetic rigidity settings of the spectrometer (in 
the range 1.3-2.0 Tm), summation over all ionic charge 
states (with corrections applied for missing charge states 
[24j]), and, finally, normalization for beam current and 
target thickness, provided fragment yield distributions 
with respect to Z, A and velocity. Further details of the 
analysis procedure can be found in |25l and in previous 
work with heavier beams ptiL I2H |28|. The yield distri- 
butions, summed over velocities, were used to obtain the 
yield ratios R 2 i{N, Z) = Y 2 {N, Z)/Y 1 (N, Z) employed in 
the present isotopic scaling studies. 



III. RESULTS AND DISCUSSION 

Before embarking on the discussion of the yield ratios 
and scaling, we will first examine the velocity and N/Z 
characteristics of the reaction products. It is well es- 
tablished that reactions between massive nuclei around 
the Fermi energy [l8| proceed via a deep inelastic trans- 
fer mechanism involving substantial nucleon exchange 
[29T I30L I3TL l32| . This mechanism is responsible for the 
creation of highly excited primary fragments that deex- 
cite to produce the observed fragments. To obtain infor- 
mation about the excitation energy of the primary frag- 
ments from the present reactions, we will first examine 
the correlation of the measured velocity with the atomic 
number. Fig. 1 presents the average velocities of the 
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fragments as a function of Z. Closed symbols correspond 
to the reactions with the neutron-rich targets and open 
symbols to those with the neutron-poor targets. In this 
figure, we observe that for fragments close to the projec- 
tile, the velocities are slightly below that of the projectile, 
corresponding to very peripheral, low-excitation energy 
events. A monotonic decrease of velocity with decreas- 
ing Z is observed, indicative of lower impact parameters, 
higher momentum transfers, and thus higher excitation 
energies. 
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FIG. 1: Average velocity versus atomic number Z 
correlations for projectile residues from the reactions of 
(a) 86 Kr(25MeV/nucleon) with 124 Sn and 112 Sn, and (b) 
86 Kr(25MeV/nucleon) with 64 Ni and 58 Ni. Full circles rep- 
resent the data with the neutron-rich targets 124 Sn and 64 Ni, 
and open circles those with the neutron-poor targets 112 Sn 
and 58 Ni. The dashed line (marked "PR") gives the velocity 
of the projectile, whereas the arrows indicate the minimum 
average residue velocities observed (see text). 

For the 86 Kr + 124 < 112 Sn reactions (Fig. la) , the de- 
scending velocity-Z correlation continues down to Z^28; 
for lower Z's, the velocity appears to increase with de- 
creasing Z. The appearance of a minimum velocity for 
Z~28 can be understood by assuming that these residues 
originate from primary fragments with a maximum ob- 
served excitation energy. Fragments with Z near the 
projectile down to Z^28 originate from evaporative type 



of deexcitation which does not modify, on average, the 
emission direction of the residues. Thus, the residue 
velocities can provide information on the excitation en- 
ergy. Residues with lower Z arise, as we will discuss 
below, from primary fragments undergoing cluster emis- 
sion and/or multifragmentation and the velocity of the 
inclusively measured fragments is not monotonically cor- 
related with the excitation energy, the mass A or the 
atomic number Z. For the 86 Kr + 64 - 58 Ni reactions, a 
similar behavior is observed. However, the decreasing 
velocity-Z correlation is observed only down to Z~32. 
It should be pointed out that fragments from this re- 
action were measured mostly outside the grazing angle, 
so that they correspond to more damped collisions, in 
such a way that the final residues receive a larger re- 
coil during the deexcitation stage and appear within this 
angular range. For Z~30-32, we observe a minimum ve- 
locity and for lower Z's an increase of the velocity with 
decreasing Z. The average velocities from the reactions 
with the neutron-rich targets are, within the experimen- 
tal uncertainties, almost the same as the corresponding 
from the reactions with the neutron-deficient targets for 
both pairs of systems. For both reactions, the ascending 
part of the velocity vs Z correlation for the lower part 
of the Z range is primarily due to the combined effect 
of angle and magnetic rigidity selection. The forward 
angle range (2.7°-5.4°) selects either the forward or the 
backward kinematical solution in the moving frame of the 
quasiprojectile undergoing cluster-emission or multifrag- 
mentation, whereas the magnetic rigidity range (1.3-2.0 
T m) subsequently selects the forward solution. 

Employing the observed minimum velocities for the 
Kr+Sn and Kr+Ni reactions and, furthermore, apply- 
ing two-body kinematics and equal division of excita- 
tion energy (which is a reasonable assumption for nearly 
symmetric systems at this energy regime |32l |33|), we 
can estimate an average excitation energy per nucleon 
for the hot quasiprojectile fragments of E*/A=2.2±0.1 
MeV for both Kr+Sn and Kr+Ni systems. Using the 
Fermi gas relationship E* = ^T 2 , with T the tempre- 
tature and K the inverse level density parameter, taken 
as K=13 MeV |34j], we can estimate the temperature of 
Kr-like quasiprojectiles as T=5.3±0.2 MeV for both sys- 
tems. This average temperature is in the range typical for 
reaction processes near the multifragmentation threshold 
[Tfll34l|. In addition, this temperature is rather close to 
the limiting temperature of excited nuclei of A~90 ac- 
cording to the recent systematics of Natowitz et al. [35| . 
Since these excitation energy estimates are average val- 
ues obtained from average residue velocities, we may as- 
sume that surviving residues with even lower velocities 
originating from even higher excitation energies are pro- 
duced in the present measurements. However, due to the 
inclusive nature of the measurements and the velocity 
fluctuations due to evaporation, we cannot select residue 
source velocities lower than those indicated by the aver- 
age velocities. 

In Fig. 2, we present the average N/Z values for each 
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Z for the Kr+Sn reactions (upper panel) and Kr+Ni 
reactions (lower panel). In the figures, the N/Z val- 
ues from the reactions involving the neutron-rich tar- 
gets are shown with full circles and those involving the 
neutron-poor targets with open circles. The horizontal 
dashed line (marked "PR") represents the N/Z of the 
86 Kr projectile. The thin solid line gives the line of 
f3 stability (marked "SL") obtained from the relation: 
Z/3 = A/ (1.98 + 0.0155A 2 / 3 ) [36|. Finally, the dotted 
line represents the evaporation attractor line (EAL) j3?J 
corresponding to the locus of fragment yields produced 
by evaporation of neutrons and charged particles from 
highly excited nuclei close to stability. (The statistical 
deexcitation code GEMINI was used to obtain the posi- 
tion of the EAL [H^O 
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FIG. 2: Average N/Z versus atomic number Z cor- 
relations for projectile residues from the reactions of 
(a) 86 Kr(25MeV/nucleon) with 124 Sn and 112 Sn, and (b) 
86 Kr(25MeV/nucleon) with 64 Ni and 58 Ni. Full circles rep- 
resent the data with the neutron-rich targets 124 Sn and 64 Ni, 
and open circles those with the neutron-poor targets 112 Sn 
and 58 Ni. The dashed line (marked "PR") gives the N/Z of 
the projectile, the solid line (marked "SL") is the line of sta- 
bility and the dotted line (marked "EAL" ) is the evaporation 
attractor line p37| (see text). 

For the Kr+Sn reactions, we observe that the average 
N/Z of the fragments close to the projectile rapidly de- 



creases with deceasing Z down to approximately Z=26- 
28. For lower Z, the average N/Z appears to be grossly 
constant with a possible dip around Z=15. The Z value 
below which the average N/Z of the fragments is roughly 
constant coincides with the Z value corresponding to the 
minimum residue velocity, as previously discussed in re- 
lation to Fig. 1. As expected, the N/Z of the frag- 
ments from the reactions with the neutron-rich 124 Sn are 
larger than those from the reactions with the neutron- 
poor 112 Sn target. This difference appears to be largest 
below Z~26, namely close to the onset of cluster emission 
and/or multifragmention. In this region, the fragment 
N/Z deviates from its course towards or beyond (to the 
neutron-poor side of) the stability line and appears to 
be roughly constant. This feature of multifragmentation 
of very neutron-rich heavy systems may be exploited in 
the production of neutron-rich rare isotopes and it will 
be the topic of a subsequent study. 

For the Kr+Ni reactions, similar observations can be 
made in Fig. 2b. For these reactions, the rapid decrease 
of N/Z is observed down to Z=32, corresponding to the 
Z with the minimum average velocity (Fig. lb). Below 
Z=32, the N/Z shows a slight decrease and then, below 
Z=26 it becomes roughly constant (again with a dip at 
Z~15, as for Kr+Sn). 

Having examined the excitation energy and N/Z char- 
acteristics of the measured residue data, we turn our dis- 
cussion to the scaling properties of the fragment yields. 
According to the equilibrium limit of the grand-canonical 
ensemble, a thermally equilibrated system undergoing 
statistical decay can be characterized by a primary frag- 
ment yield with neutron number N, and proton number 
Z of the form HHIl: 



Y(N, Z) = F(N, Z) exp[B{N, Z)/T] exp(Nfx n /T+Z(j, p /T) 

(1) 

where the factor F(N, Z) represents contributions due 
to the secondary decay from particle stable and unstable 
states to the final ground state; /x ra and \x v are the neutron 
and proton chemical potentials; B(N, Z) is the ground 
state binding energy of the corresponding fragment, and 
T is the temperature. 

A direct comparison of the observed yield distributions 
with Eq. 1 is not possible due to the distortions intro- 
duced by the sequential decay of the primary fragments. 
However, it has been shown that the ratio of the 
yields Y 2 (N, Z) /Y X (N, Z) of a given fragment (N,Z) from 
two different reacting systems, having similar masses and 
excitation energies, but differing only in N/Z can be used 
to obtain information about the excited primary frag- 
ments. If the main difference between the two systems 2 
and 1 is the isospin, then the binding energy terms in Eq. 
1 cancel out in the ratio Y 2 (N, Z)/Y 1 (N, Z). Further- 
more, if one assumes that the influence of the secondary 
decay on the yields is similar for the two reactions, then 
a relation of the form 
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R 21 {N, Z) = y 2 (JV, Z)/Y 1 {N, Z) = Cexp(aN + (3Z) 

(2) 

can be obtained with a = A/i„/T and /3 = A/z p /T, 
with A//„ and A^ p being the differences in the neutron 
and the proton chemical potentials of the fragmenting 
systems. C is an overall normalization constant. The 
ratio R,2i(N, Z) is insensitive to the sequential decay and 
thus it can provide information on the decaying excited 
primary fragments. 

From the present data, we construct the yield ratio 
R.2i(N,Z) using the convention that index 2 refers to 
the more neutron-rich system and index 1 to the less 
neutron-rich one. Fig. 3 shows the yield ratios R2i(N,Z) 
as a function of fragment neutron number N for selected 
isotopes (top panel) and proton number Z for selected 
isotones (bottom panel) for the Kr+Sn reactions. The 
different isotopes and isotones considered are shown by 
alternating filled and open symbols for clarity. In Fig. 
4, the corresponding ratios for the Kr+Ni reactions are 
shown. 

As it can be seen in the top panels of Figs. 3 and 4, the 
ratios for each element Z exhibit a remarkable exponen- 
tial behavior. For each element, an exponential function 
of the form Cexp(a N) was fitted to the data and also 
shown in Fig. 4 for the selected elements. In the semi-log 
representation, the lines for each element are nearly par- 
allel up to Z-28 for Kr+Sn and up to Z-34 for Kr+Ni 
(in the latter case exhausting the whole range of observed 
fragments). For heavier fragments from the Kr+Sn sys- 
tems, the fits to the data show gradual decrease in the 
slopes with increasing Z of the fragments. An analogous 
behavior is observed in the lower panels of Figs. 3 and 
4, where the ratio R,2i(N,Z) is plotted as a function of 
Z for various isotones. The solid lines correspond to ex- 
ponential fits using the expression Gexp(pZ) . The ratios 
R,2i(N,Z) appear to lie along straight lines with nearly 
similar negative slopes (3, for all isotones up to N~34 for 
Kr+Sn and N~44 for Kr+Ni (again, in the latter case 
covering the whole fragment range). For larger N from 
Kr+Sn, the slopes f3 increase with increasing N of the 
fragments. The positive slopes in the upper panels of 
Figs. 3 and 4 indicate that neutron-rich fragments are 
more efficiently produced, as expected, from the more 
neutron-rich systems. Similarly, the negative slopes in 
the bottom panels of these figures indicate that proton- 
rich fragments are more efficiently produced from the 
more proton-rich systems. 

In Fig. 5, we present the slope parameters a (upper 
panel) and (3 (lower panel) of the exponential fits (ob- 
tained as described for Figs. 3 and 4) as a function of 
Z and N. For the Kr+Sn reactions, the slope parameter 
a is roughly constant with an average value of 0.43 in 
the range Z=10-26 and then decreases for Z>26. Simi- 
larly, the parameter (3 appears to be nearly constant at 
an average value of -0.50 up to N~34 and rise for N>34. 
Consequently, for the Kr+Sn systems under the present 
experimental conditions, isoscaling (as defined in 
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FIG. 3: Yield ratios R 2 i(N, Z) = Y 2 (N, Z)/Yi{N, Z) of pro- 
jectile residues from the reactions of 86 Kr(25MeV/nucleon) 
with 124,112 Sn (a) with respect to N for the Z's indicated, and 
(b) with respect to Z for the N's indicated. The data are 
given by alternating filled and open circles, whereas the lines 
are exponential fits (see text). 



holds up to Z~26 and N~34. This fragment range cor- 
responds to primary events with the maximum observed 
excitation energy of 2.2 MeV/nucleon and temperature 
of 5.3 MeV, as estimated previously. For the Kr+Ni re- 
actions, the slope parameter a appears to be constant in 
the whole range Z=10-36 at an average value a=0.27. 
Similarly, the parameter (3 remains constant at an aver- 
age value of -0.34 in the whole N range. In the case of 
the Kr+Ni systems under the present experimental con- 
ditions (observation outside the grazing angle), isoscaling 
holds in the whole range of observed fragments, corre- 
sponding to primary events having on average excitation 
energy of 2.2 MeV/nucleon and temperature of 5.3 MeV. 
The near constancy of the isoscaling parameters a and (3 
(except for near-projectile fragments for the Kr+Sn sys- 
tem) corroborates the statistical nature of the fragment 
production under the excitation energy and temperature 
values extracted from the minimum observed residue ve- 
locities. In general, the extracted isoscaling parameters 
a and (3 are in agreement with the corresponding values 
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FIG. 4: Yield ratios R 2 i(N, Z) = Y 2 (N, Z)/Y X (N, Z) of pro- 
jectile residues from the reactions of 86 Kr(25MeV/nucleon) 
with 64,58 Ni (a) with respect to N for the Z's indicated, and 
(b) with respect to Z for the N's indicated. The data are 
given by alternating filled and open circles, whereas the lines 
are exponential fits (see text). 



reported in the literature (e.g. 0,0). 

In the following discussion, we will use the isoscaling 
parameters herein obtained to extract information about 
the hot primary fragments. As we have already men- 
tioned, the parameters a and /3 are related to the dif- 
ferences in the neutron and proton chemical potentials 
of the corresponding primary fragmenting systems via 
the relations: a = A/i n /T and (3 — A/i p /T. Using the 
average values of the parameters a and f3 and the esti- 
mated temperatures, we can obtain the differences in the 
chemical potentials, as summarized in Table I. As seen 
in the table, these values are close (apart from a differ- 
ence in sign) 0, 0] to the corresponding differences in 
the neutron and proton separation energies, as obtained 
from mass tables ^lj. The corresponding values for the 
protons appear to be remarkably close (Table I). 

It has been shown 0,^1 that the isoscaling parameter 
a is directly related to the coefficient C sym of the sym- 
metry energy term of the nuclear binding energy. The 
following relation has been obtained both in the frame- 
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FIG. 5: (a) Isoscaling parameter a as a func- 

tion of Z for projectile residues from the reactions 
86 Kr(25MeV/nucleon)+ 124 ' 112 Sn (closed circles) and 
86 Kr(25MeV/nucleon)+ 64 ' 58 Ni (open circles). The 
straight lines are constant value fits for each system 
(see text), (b) Isocaling parameter j3 as a function of N 
for 86 Kr(25MeV/nucleon)+ 124 ' 112 Sn (closed circles) and 
86 Kr(25MeV/nucleon)+ 64,58 Ni (open circles). The straight 
lines are again constant value fits for each system. 



work of the grand canonical limit of the statistical multi- 
fragmentation model (SMM) [^| and in the expanding- 
emitting source (EES) model [14): 



a = 4 ■ 



a 



syrn 



((^) 2 -(^) 2 ) 

Ax Ai 



(3) 



where Zi,A\ and ^2,^2 refer to the fragmenting 
quasiprojectiles from reactions 1 and 2 respectively. Us- 
ing this relation, the extracted values of a and T for the 
fragments corresponding to multifragmentation events 
and assuming N/Z equilibration (see following discussion 
and Fig. 6), we get for C sym the values of 27.2 for Kr+Sn 
and 23.1 for Kr+Ni (Table I). These estimates of the 
symmetry-energy coefficient are in reasonable agreement 
with the standard value C sym = 25 MeV [l5| . 

In principle, Eq. [3] can serve as the basis for deter- 
mining C sym for expanded multifragmenting nuclei & 
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TABLE I: Excitation energy E*/A (MeV/nucleon), temper- 
ature T (MeV), isoscaling parameters a and /3, differences 
in neutron and proton chemical potentials A/i n (MeV) and 
A/ip (MeV), differences in neutron and proton separation 
energies AS n (MeV) and AS P (MeV), and extracted sym- 
metry energy coefficient C sym (MeV) for the systems 86 Kr 
(25 MeV/nucleon) + 124 ' 112 Sn, and 86 Kr (25 MeV/nucleon) 
_l_ 64,58-^ ^ gee ^ext). Parentheses give estimated errors (one 
standard deviation). 



Rxn E*/A T 



(3 A/in A/Hp ASn ASp Csyn 



Kr+Sn 2.2 5.3 0.43 -0.51 2.28 -2.65 -1.44 2.23 27.2 
(0.1) (0.2) (0.03) (0.03) (0.18) (0.19) (2.2) 




15 20 25 30 
Atomic Number Z 



Kr+Ni 2.2 5.3 0.27 -0.34 1.43 -1.80 
(0.1) (0.2) (0.02) (0.02) (0.11) (0.13) 



FIG. 6: Difference in N/Z of the fragmenting quasiprojectiles 
2.75 1.77 23.1 (obtained using Eq. ^ as a function of Z for projectile residues 
(1-9) from the reactions 86 Kr(25MeV/nucleon)+ 124,112 Sn (closed 
circles). The dashed line gives the N/Z difference of isospin- 
equilibrated fragmenting quasiprojectiles. 



thus probing the density dependence of the symmetry 
energy term of the nuclear equation of state. However, 
such experimental determination requires systematic iso- 
topically resolved yield measurements with good exci- 
tation energy characterization. With the present spec- 
trometric approach, in order to avoid the uncertainties 
in excitation energy division, as well as in N/Z equili- 
bration, such measurements may be performed for sym- 
metric pairs differing in isospin (e.g. 64 Ni+ 64 Ni and 

58 M+ 58 M Qr 124 Sn+ 124 Sn &nd U2 Sn+ U2 Sn )_ J n such 

cases, the excitation energies (and temperatures) can be 
obtained from the residue velocities and the average N/Z 
of the quasiprojectiles can be taken as that of the reacting 
system. 

Finally, we will employ Eq. E|to provide an interpreta- 
tion of the variation of the isoscaling parameter a with Z 
for the Kr+Sn reactions and relate this variation to the 
isospin equilibration process. After some manipulation, 
from Eq. [3] we obtain: 

« = 8^(§)L e A(§V (4) 

where (Z/A) ave is the average Z/A of the quasipro- 
jectiles, taken to be the average Z/A of the composite 
systems 86 Kr+ 124 Sn and 86 Kr+ 112 Sn, and A(N/Z) qp ex- 
presses the N/Z difference of fragmenting quasiprojectiles 
corresponding to a given value of fragment Z, isoscaling 
parameter a and temperature T. Assuming that frag- 
ments are produced at normal density, using C sym = 25 
MeV [13, the a values obtained from the isoscaling fits 
and the temperatures determined from excitation ener- 
gies (obtained from residue velocities) , we can determine 
the value of A(N/Z) qp as a function of observed fragment 



Z, as shown in Fig. 6. In this figure, the horizontal line 
expresses the N/Z difference of fragmenting quasiprojec- 
tiles under the condition of isospin equilibration. 

The essentially monotonic increase of A(N/Z) qp with 
decreasing Z can be understood as follows. Fragments 
close to the projectile are produced in very peripheral col- 
lisions in which a small number of nucleons is exchanged 
and thus, the N/Z difference of the fragmenting quasipro- 
jectiles from 86 Kr+ 124 Sn and 86 Kr+ 112 Sn is small. Frag- 
ments progressively further from the projectile originate 
from collisions with larger projectile-target overlap in 
which a large number of nucleons is exchanged. In such 
cases, the N/Z difference of the fragmenting quasipro- 
jectiles is progressively larger, eventually approaching, 
for the present energy regime, the N/Z difference corre- 
sponding to isospin equilibration. Consequently, the vari- 
ation of the isoscaling parameter a with Z (and, equiva- 
lently, (3 with N) for near-projectile residues is a result of 
the variation in the N/Z of the fragmenting quasiprojec- 
tiles, directly related to the course towards isospin equi- 
libration. As such, the parameter a along with tempera- 
ture (excitation energy) determination can be used as a 
sensitive probe to study the isospin equilibration process 
in collisions of isospin asymmetric massive nuclei at in- 
termediate energies. A detailed account on this approach 
is reported elsewhere [42^. 

The present high-resolution mass spectometric data 
provided information on the isoscaling properties of 
heavy projectile fragments and, subsequently, informa- 
tion on the de-excitation of hot primary fragments, as 
well as, the degree of N/Z equilibration. Along these 
lines, it should be noted that with the same spectromet- 
ry techniques, we recently obtained information on the 
details of the de-excitation mechanisms of heavy frag- 
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ments and the onset of multifragmentation in the reac- 
tion of 124 Sn(25MeV/nucleon) with 27 A1 43]. These re- 
cent efforts in the Fermi energ y reg ime (following earlier 
work with heavier beams |2(i l27j 1. along with similar 
studies at relativistic energies (e.g. HH an d refer- 
ences therein) indicate that detailed isotopic and kine- 
matic investigation of heavy residues provides essential 
information on the reaction mechanisms and the details 
of the de-excitation processes that can complement the 
results of full-acceptance devices. As a possible next step, 
combination of high-resolution mass separators with ad- 
vanced large-acceptance multidctcctor arrays may allow 
novel and exciting studies of reaction dynamics, further- 
more taking advantage of the availability of rare isotope 
beams. 



IV. SUMMARY AND CONCLUSIONS 

In summary, the isotopic scaling of heavy projectile 
residues from the interaction of 25 MeV/nucleon 86 Kr 
projectiles with 124 Sn, 112 Sn and 64 Ni, 58 Ni targets has 
been studied. Isotopically resolved yield distributions of 
projectile fragments in the range Z= 10-36 from the above 
reaction pairs were measured with the MARS recoil sep- 
arator in the angular range 2.7°-5.4°. For these deep in- 
elastic collisions, the velocities of the residues, monoton- 
ically decreasing with Z down to Z~26-28, are employed 
to characterize the excitation energy. The yield ratios 
R,2i(N,Z) for each pair of systems are found to exhibit 
isoscaling, namely, an exponential dependence on the 
fragment atomic number Z and neutron number N. The 
isoscaling is found to occur in the residue Z range corre- 
sponding to the maximum observed excitation energies. 
The corresponding isoscaling parameters are a=0.43 and 



(3— -0.50 for the Kr+Sn system and a— 0.27 and (3=- 
0.34 for the Kr+Ni system. For the Kr+Sn system, for 
which the experimental angular acceptance range lies in- 
side the grazing angle, isoscaling was found to occur for 
Z<26 and N<34. The values of the isoscaling parame- 
ters are used to extract the neutron and proton chemical 
potentials, as well as obtain an estimate of the symmetry 
energy coefficient of highly excited primary fragments. 
For heavy fragments from Kr+Sn, the isoscaling behav- 
ior breaks down and the parameters are found to vary 
monotonically, a decreasing with Z and (3 increasing with 
N. This variation of the isoscaling parameters for near- 
projectile residues is found to be related to the evolution 
towards isospin equilibration and, as such can serve as 
a tracer of N/Z equilibration. The present isotopically 
resolved heavy-residue data demonstrate the occurence 
of isotopic scaling from the intermediate mass fragment 
region to the heavy-residue region. Such high-resolution 
mass spectrometric data have the potential to provide 
important information on the role of isospin and isospin 
equilibration in peripheral and mid-peripheral collisions, 
complementary to that accessible from advanced large- 
acceptance multidetector devices. 
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